Electro-Osmotic Flow and Mass Transfer through a Rough Microchannel with a Modulated Charged Surface

In this paper, we investigate the electro-osmotic flow (EOF) and mass transfer of a Newtonian fluid propelled by a pressure gradient and alternating current (AC) electric field in a parallel microchannel with sinusoidal roughness and modulated charged surfaces. The two-wall roughness is described by in-phase or out-of-phase sine functions with a small amplitude δ. By employing the method of perturbation expansion, the semi-analytical solutions of the Poisson–Boltzmann (P–B) equation based on the Debye–Hückel approximation and the modified Navier–Stokes (N–S) equation are obtained. The numerical solution of the concentration equation is obtained by the finite difference method. The effects of sinusoidal roughness, modulated charged surface, and the AC electric field on the potential field, velocity field, and concentration field are discussed. Under the influence of the modulated charged surface and sinusoidal roughness, vortices are generated. The velocity oscillates due to the effect of the AC electric field. The results indicate that solute diffusion becomes enhanced when the oscillation Reynolds number is below a specific critical value, and it slows down when the oscillation Reynolds number exceeds this critical value.

Electro-osmosis has become the preferred method for fabricating microdevices due to its low energy consumption [19,20].When the electrolyte solution is filled in charged surface microchannels, ions with opposite charges to the wall are attracted to the charged surface, thus forming an electric double layer (EDL) [13,17].When the electric field is applied to both ends of the microchannel, the ions in the electrical double layer (EDL) begin to move under the influence of the electric field force.Because the fluid is viscous, it flows along with the ions in the microchannel.This flow is called the electro-osmotic flow (EOF) [6].Chang et al. [21] investigated the electro-osmotic flow in a microchannel with a charge slip on a plate surface.The influence of the thermodiffusive effect on the local Debye length thickness in a purely electro-osmotic flow in a parallel flat plate microchannel was theoretically studied by Hernández et al. [22].
Previous studies have been based on uniform potential.However, due to the limitation of the industrial level, it is not possible to achieve uniform distribution of surface charges [23,24].Mandal et al. [25] used a modulated charged surface to simulate charge inhomogeneity to study the electro-osmotic flow of superimposed fluids under narrow constraints.Bian and Li [26] consider fluid flow in a microchannel with a modulated charged surface.Akhtar Shehraz et al. [27] used fractional derivatives to study the electro-osmotic flow of Maxwell fluids in microchannels with asymmetric wall potential.Wang and Li [28] studied the EOF and heat transfer in the polyelectrolyte-grafted microchannels with modulated charged surfaces.The above results show that the modulation of charged surfaces can generate eddy currents in the flow field, which have positive effects on fluid mixing, solute diffusion, and heat transfer.
Given that the current manufacturing process is unable to achieve perfect smoothness, there will always be some degree of wall roughness.Sadia Siddiqa [29] conducted a study on the transient analysis of natural convective heat transfer based on vertical wavefronts and obtained the solution to the equation using the coordinate transformation method.Buren et al. [30] analyzed the critical wave value of the influence of wall roughness on velocity and potential distribution using the perturbation expansion method.Li et al. [31] investigated the electromagnetohydrodynamic (EMHD) flow in a three-dimensional corrugated wall microchannel.Chang et al. [32] discussed the impact of sinusoidal roughness on the AC EOF of Maxwell fluids in parallel microchannels.
The diffusion process of solutes is the mass transport due to the inhomogeneity of molecular diffusion and fluid flow velocity.The problem of mass transport in microfluidic systems has important guiding significance in practical applications such as drug delivery or toxin separation in medical care and biological systems [33].More recently, the applications of the AC electric fields have been extended to electro-osmotic flow with enhanced mass transfer, mixing, and material separation [34,35].Medina et al. [36] studied the pulsating electro-osmotic flow and solute diffusion in microchannels with wall potential asymmetry.Li and Jian [37] studied the problem of solute diffusion in polyelectrolytegrafted nanochannels driven by the AC electric field.They discovered that there is a critical oscillation Reynolds number value for solute diffusion driven by the AC electric field.The electromagnetic hydrodynamic flow and mass transfer in curved rectangular microchannels were studied by Liu and Jian [38].They found that the effective diffusivity increases with the oscillation Reynolds number.
According to the author's knowledge, currently, no relevant studies have discussed the effects of roughness and modulated charged surface coupling on fluid flow.Therefore, this paper examines the electro-osmotic flow in a microchannel with a rough and modulated charged surface, driven by both the AC electric field and the pressure gradient.The P-B equation based on the Debye-Hückel approximation, and the improved N-S equation are solved using the asymptotic expansion method.The numerical solution of the concentration equation using perturbation expansion is obtained by the finite difference method.Finally, the effects of sinusoidal roughness, modulated charged surface, and the AC electric field on potential, velocity, and mass transport are discussed in detail.

Mathematics Model
In this study, we consider the alternating current electro-osmotic flow and mass transfer of Newtonian fluid in a microchannel with sinusoidal roughness and modulated charged surfaces.Shown as in Figure 1, the length of the mathematical model in the x * -axis direction is denoted as L (L ≫ 2H).First, let us consider the low Reynolds number flows in the microchannels.The electric field E x * (t * ) is aligned with the x * -axis.The fluid is propelled by the AC electric field and pressure gradient.The positions of the upper and lower sinusoidal roughness walls (y * u and y * l ) can be described by periodic sine waves as follows: where λ* is the wave number, and δ ≪ 1 is the amplitude.y * + l = −H + δH sin(λ * x * ) denotes the in-phase walls, which means that the lower and upper walls are not symmetric.y * − l = −H − δH sin(λ * x * ) denotes the out-of-phase walls, which means that the lower and upper walls are symmetric.
where λ* is the wave number, and δ ≪ 1 is the amplitude.y l * + = -H + δHsin λ * x * denotes the in-phase walls, which means that the lower and upper walls are not symmetric.y l *-= -H -δHsin λ * x * denotes the out-of-phase walls, which means that the lower and upper walls are symmetric.

Velocity Distribution
The microchannel walls adopt an asymmetric modulated potential, and the potential distribution function ψ * x * , y * of the upper and lower walls is expressed as follows: where ξ u * and ξ l * are the amplitudes of the upper and lower walls, and α and β are constants.It is assumed that the electrolyte solution fills the microchannel.The potential ψ * x * ,y * and the net charge density  * can be described by the P-B equation [37]: where  * = 2 zesinh(zeψ * / k B T a ) is Bol mann distribution.The symbol  represents the concentration of ions in the solution.z represents the valence of the ion.e is the charge of the electron, T a stands for the absolute temperture, k B denotes the Bol mann constant, and ε signifies the dielectric constant.By the Debye-Hückel approximation, we obtain the term sinh(zeψ * / k B T a ) ≈ zeψ * / k B T a .From this we can obtain the approximate P-B equation: where κ = 1/[(ek B T a )/(2   )] 1/2 denotes the thickness of the EDL, also known as the Debye length.The boundary conditions are (2a) and (2b).Since both sinusoidal roughness and modulating charge surface effects are taken into account, the velocity is two-dimensional, i.e.,  * = (U * ,  * , 0).Then, the continuity equation and the N-S equation can be expressed as [36,37]:

Velocity Distribution
The microchannel walls adopt an asymmetric modulated potential, and the potential distribution function ψ * (x * , y * ) of the upper and lower walls is expressed as follows: where ξ * u and ξ * l are the amplitudes of the upper and lower walls, and α and β are constants.It is assumed that the electrolyte solution fills the microchannel.The potential ψ * (x * , y * ) and the net charge density ρ * e can be described by the P-B equation [37]: where ρ * e = 2n 0 zesinh(zeψ * / k B T a ) is Boltzmann distribution.The symbol n 0 represents the concentration of ions in the solution.z represents the valence of the ion.e is the charge of the electron, T a stands for the absolute temperture, k B denotes the Boltzmann constant, and ε signifies the dielectric constant.By the Debye-Hückel approximation, we obtain the term sinh(zeψ * / k B T a ) ≈ zeψ * / k B T a .From this we can obtain the approximate P-B equation: where 1/2 denotes the thickness of the EDL, also known as the Debye length.The boundary conditions are (2a) and (2b).Since both sinusoidal roughness and modulating charge surface effects are taken into account, the velocity is two-dimensional, i.e., U * = (U * , V * , 0).Then, the continuity equation and the N-S equation can be expressed as [36,37]: There is an external alternating electric field E x * (t) along the x * -axis, and the external force acting on the fluid is F = ρ * e E x * (t * ).Suppose that the AC electric field can be expressed in complex form: where Re{} represents the real part of the complex number, E 0 is the amplitude of the applied AC electric field, i is the imaginary unit, and ω is the oscillating angular frequency of the applied AC electric field.The momentum Equations ( 5) and ( 6) are written in a two-dimensional component form.The final Equations ( 5) and ( 6) are written as [31]: The boundary conditions for Equations ( 8)-( 10) are as follows: Based on the expression of the AC electric field, we assume the complex form of the velocity [37]: where u * (x * , y * ) and v * (x * , y * ) are the amplitudes of the periodic velocities in the x * and y * directions, respectively.The following non-dimensional variables are introduced: where UHS is the Helmholtz velocity of the Newtonian fluid.The P-B equation, the continuity equation, and the N-S equation can be transformed as: ∂u ∂x − ∂P ∂y where Re ω = ρωH 2 µ denotes the oscillation Reynolds number and the competition between the characteristic diffusion time and the characteristic time associated with the oscillatory electric field [37].The non-dimensional boundary conditions for Equations ( 14)- (17) are u(x, y) ± = 0 at y = y u and y = y ± l , (18c) v(x, y) ± = 0 at y = y u and y = y ± l , where boundary conditions (18c,18d) illustrate the no-slip and non-permeable boundary conditions.The asymptotic method applies to problems in which the geometry of the microchannel is complicated and irregular in the sense that the boundaries do not correspond to coordinate surfaces of any known analytic coordinate system, but are nevertheless "near" to such coordinate surfaces [39].For small values of δ ≪ 1, function (could be u, v, ψ, P, φ, C) can be expanded asymptotically by means of perturbation [30,31]: Equation ( 19) is expanded up to the second-order terms.For the asymptotic expansion of the boundary condition (18a-18d), the Taylor expansion at the mean wall position y = 1 and y = −1 yields the following: The asymptotic expansion equation of (u, v, ψ, P) is substituted into Equations ( 14)-( 17) and (20a-20f) to obtain the zero-order, first-order, and second-order equations, as well as the boundary conditions.

The Zeroth-Order Equations
The zeroth-order equations are as follows: The zeroth-order boundary conditions are as follows: v 0 (x, y) = 0 at y = 1 and y = −1. (25d) The Equation ( 21) and boundary conditions (25a,25b) are used to solve the zerothorder potential solution.Based on the superposition principle, the solution of the potential is as follows: The expression of coefficient is in Appendix A. In the following calculation, in order to solve the two-dimensional velocity, we introduce the zeroth-order stream function φ 0 according to the zeroth-order continuity equation: By substituting Equation ( 27) into Equations ( 22)-( 24) of zero-order, the pressure P 0 is eliminated, and the equation of φ 0 is obtained as follows: The solution of the zeroth-order stream function has the form as follows: The expression of coefficient is in Appendix A. Then, the zeroth-order velocity u 0 , v 0 can be obtained through Equation (31).

The First-Order Equations
The first-order equations are obtained as follows: The first-order boundary conditions are as follows: According to the superposition principle and boundary condition (34a), the first-order potential solution has the following form: The expression of coefficient is in Appendix A. Introducing first-order stream function φ 1 into Equations ( 31)-( 33), the pressure P 1 can be eliminated, yielding the following equation: The analytical solution of the first-order stream function φ 1 is in the following form: The expression of coefficient is in Appendix A. Then, the first-order velocity u 1 , v 1 can be obtained.

The Second-Order Equations
As in the above process, the second-order equations are obtained as follows: The second-order boundary conditions are as follows: According to the superposition principle and considering Equations (38) and the boundary conditions (42a,42b), the second-order potential solution has the following form: The expression of coefficient is in Appendix A. Introducing first-order stream function φ 2 into Equations ( 39)-( 41), the pressure P 2 can be eliminated, yielding the following equation: The analytical solution of the first-order stream function φ 1 is in the following form: The expression of coefficient is in Appendix A. Then, the second-order velocity u 2 , v 2 can be obtained.
From Sections 2.1.1-2.1.3,we can obtain the semi-analytical solution of the velocity of the form as follows: (46)

Concentration Distribution
In this study, we assume that a container with a constant concentration C * 0 of substance is placed at one end of the microchannel.As stated in Section 3.1, the velocity is twodimensional; therefore, solute diffusion in a microchannel can be described by a twodimensional unsteady convection-diffusion equation [36][37][38]: where C * (x * , y * , t * ) is the concentration and D is the diffusion.The boundary conditions are as follows: Equation (48a) indicates that the concentration at the initial position of the microchannel is a constant C * 0 at any time.Equation (48b) indicates that the concentration in the channel is zero at the initial time.Equation (48c) is the non-permeable boundary condition.The unsteady convection-diffusion equation and the boundary conditions are transformed into nondimensional forms by employing H/U HS and C * 0 as time and concentration scale ( t = t * t 0 , t 0 = H U HS , C = C * C 0 , respectively.The dimensionless concentration C(x, y, t) is governed by the following: where Pe= HU HS D denotes the Péclet number, and its physical meaning is the ratio of the rate of convection to the rate of diffusion with the following nondimensional boundary conditions: The dimensionless variable C can be expanded asymptotically by utilizing perturbation: From Equations ( 49)-( 51), the zeroth, first, and second order equations, as well as the boundary condition, are obtained.
where U 0 (x, y, t) = Re u 0 (x, y)e iωt , V 0 (x, y, t) = Re v 0 (x, y)e iωt . where where U 2 (x, y, t) = Re u 2 (x, y)e iωt , V 2 (x, y, t) = Re v 2 (x, y)e iωt .The analytical solution of Equations ( 52)-( 54) is still difficult to obtain.Therefore, to obtain the numerical solution of the equation, the finite difference method (FDM) [40] is utilized in this study.The detailed finite difference scheme and algorithm can be found in Appendix B.

Results and Discussion
The diffusion of solutes at one end of the microchannel, the microchannel with the rough modulated charged surface driven by an applied electric field, is studied.Some dimensionless parameters play an important role in the fluid motion in the microchannel, such as the oscillation Reynolds number Re ω , the roughness amplitude δ, the modulation parameters α, β, and so on.Based on the previous studies, the physical parameters are as follows [32,36]: the half of the microchannel H = 100 µm, the density of the fluid ρ = 10 3 kg•m −3 , and the dynamic viscosity µ = 10 −3 kgm −1 s −1 .At the same time, the applied electric field frequency ranges from 0 to 1.6 kHz, and the variation range of the corresponding angular frequency ω is from 0 to 10 4 s −1 .Therefore, the oscillation Reynolds number Re ω could take a value between 0 and 100.In the following discussion about the upper and lower modulated charged surfaces, the ranges of the amplitudes (ξ u and ξ l ), constants (α and β), and mode frequencies (λ) are 0~0.2,0~6, and 0~4 respectively.The molecular diffusion coefficient D varies from O(10 −8 )m 2 s −1 to O(10 −9 )m 2 s −1 [36].The Péclet number Pe varies from 0.3 to 1, and the range of δ is 0 ≤ δ ≤ 0.1.The initial concentration C 0 = 10 3 mol•m −3 [41].

Electric Potential Field
Figure 2 illustrates the impact of sinusoidal roughness and modulated charged surfaces on potential fields.When λ = 0 and δ = 0 in Figure 2a, the wall represents a smooth parallel plate microchannel and the potential is uniform.At this time, the lines of potential near the wall are denser, indicating a higher potential intensity.It can also be seen from Figure 2a that the potential distribution is not uniform, which is caused by the difference between ξ u and ξ l ; that is, the two walls potentials are not symmetrical.Figure 2b,c represent potential fields on in-phase walls, and Figure 2d,e represent potential fields on out-of-phase walls.By comparing Figure 2b-e with Figure 2a, it can be observed that the potential distribution on the wall is non-uniform because of the impact of the modulated charged surface and roughness.Then, the potential distribution on the walls will cause the potential inside the microchannel to be uneven.Figure 3 shows the distribution of streamlines with the sinusoidal roughness amplitude δ, modulated parameters α, β, and Re ω under the conditions of λ = 1, κ = 2, ξ = 0.1, ξ = 0.15. Figure 3a depicts streamlines in the non-rough microchannels with asymmetric wall potentials.The flow at this point is un-directional.Figure 3b-d     When α = β = 0, δ = 0, κ = 20, ξ u = ξ l = 1, and λ = 1, Figure 4 shows a comparison between the experimental results (solid lines) of Medina et al. [36] and the semi-analytic solution (circles) presented in this paper.Four dimensionless transient velocity curves at different times are drawn in the figure, and the roughness and modulated wall potential are not considered in Figure 4.As can be seen from the figure, there is a high consistency between the two curves.When α = β = 0, δ = 0, κ = 20, ξ = ξ = 1, and λ = 1, Figure 4 shows a comparison between the experimental results (solid lines) of Medina et al. [36] and the semi-analytic solution (circles) presented in this paper.Four dimensionless transient velocity curves at different times are drawn in the figure, and the roughness and modulated wall potential are not considered in Figure 4.As can be seen from the figure, there is a high consistency between the two curves.5 illustrates that as the oscillatory Reynolds number increases, the oscillatory nature of the dimensionless velocity  ± becomes apparent.With the increase in the oscillatory Reynolds number, the amplitude of the velocity first increases and then decreases.Physically, in the case of a high oscillating frequency, the diffusion time is much longer than the oscillating period, which leads to a decrease in the velocity amplitude.Figure 6 shows that as the oscillatory Reynolds number increases, the amplitude of the vertical velocity  ± decreases.There is no significant difference in the oscillations of vertical velocity for different oscillatory Reynolds numbers.When κ = 2, ξ u =0.1, ξ l = 0.15, λ = 1, ωt = 1, δ = 0.03, and α = 3, β = 4, Figures 5 and 6 describe the variations of dimensionless velocities U ± and V ± with variable Re ω .Figure 5 illustrates that as the oscillatory Reynolds number increases, the oscillatory nature of the dimensionless velocity U ± becomes apparent.With the increase in the oscillatory Reynolds number, the amplitude of the velocity first increases and then decreases.Physically, in the case of a high oscillating frequency, the diffusion time is much longer than the oscillating period, which leads to a decrease in the velocity amplitude.Figure 6 shows that as the oscillatory Reynolds number increases, the amplitude of the vertical velocity V ± decreases.There is no significant difference in the oscillations of vertical velocity for different oscillatory Reynolds numbers.When α = β = 0, δ = 0, κ = 20, ξ = ξ = 1, and λ = 1, Figure 4 shows a comparison between the experimental results (solid lines) of Medina et al. [36] and the semi-analytic solution (circles) presented in this paper.Four dimensionless transient velocity curves at different times are drawn in the figure, and the roughness and modulated wall potential are not considered in Figure 4.As can be seen from the figure, there is a high consistency between the two curves.5 illustrates that as the oscillatory Reynolds number increases, the oscillatory nature of the dimensionless velocity  ± becomes apparent.With the increase in the oscillatory Reynolds number, the amplitude of the velocity first increases and then decreases.Physically, in the case of a high oscillating frequency, the diffusion time is much longer than the oscillating period, which leads to a decrease in the velocity amplitude.Figure 6 shows that as the oscillatory Reynolds number increases, the amplitude of the vertical velocity  ± decreases.There is no significant difference in the oscillations of vertical velocity for different oscillatory Reynolds numbers.7a,c, Re ω = 0.1, and in Figure 7b,d, Re ω = 70.With the increase in the oscillatory Reynolds number Re ω , the velocity distribution gradually exhibits oscillatory characteristics.The velocity  ± oscillates in the region close to the wall.And it tends to zero in the region far from the wall at the large oscillatory Reynolds number.This is because, at high electric field frequencies, the diffusion time is greater than the oscillation period.7a,c, Re ω = 0.1, and in Figure 7b,d, Re ω = 70.With the increase in the oscillatory Reynolds number Re ω , the velocity distribution gradually exhibits oscillatory characteristics.The velocity U ± oscillates in the region close to the wall.And it tends to zero in the region far from the wall at the large oscillatory Reynolds number.This is because, at high electric field frequencies, the diffusion time is greater than the oscillation period.
Figure 8 illustrates the effect of the charge modulation surface and the sinusoidal roughness on the velocities U ± and V ± .In Figure 8a,b, the dimensionless velocity U ± profile is presented, while Figure 8c,d displays the dimensionless velocity V ± profile.Figure 8a,b indicates that the amplitude of velocity U ± becomes enhanced with the application of modulation potential.This is because as the amplitude of the modulation potential heightens, the potential in the EDL also increases, leading to a rise in the amplitude of the velocity U ± .When the roughness increases, it can be observed from Figure 8a that the velocity U + becomes incremental.However, Figure 8b shows that the velocity U − slows down with the expansion of the roughness.In Figure 8c,d, it is observed that as the amplitude of the modulation potential increases, the amplitude of the velocity V ± also increases.The rise of the roughness highlights the oscillation characteristics of the velocity V ± .Figure 8 illustrates the effect of the charge modulation surface and the sinusoidal roughness on the velocities  ± and  ± .In Figure 8a,b, the dimensionless velocity  ± profile is presented, while Figure 8c,d displays the dimensionless velocity  ± profile.Figure 8a,b indicates that the amplitude of velocity  ± becomes enhanced with the application of modulation potential.This is because as the amplitude of modulation potential heightens, the potential in the EDL also increases, leading to a rise in the amplitude of the velocity  ± .When the roughness increases, it can be observed from Figure 8a that the velocity  becomes incremental.However, Figure 8b shows that the velocity  slows down with the expansion of the roughness.In Figure 8c,d, it is observed that as the amplitude of the modulation potential increases, the amplitude of the velocity  ± also increases.The rise of the roughness highlights the oscillation characteristics of the velocity  ± .

Concentration Field
When ξ = 0.1, ξ = 0.15, α = 3, β = 4, x = 3, λ = 1, κ = 2, t = 1, and  = 0.1, Figure 9 shows the comparison of concentrations at the fixed point x = 3 of the microchannel with three different wall surfaces: non-roughness (C 0 ), in-phase (C + ), and out-of-phase (C -).The concentration C + is higher than the concentration C 0 of the non-roughness, which indicates that the solute diffusion rate on the in-phase walls is higher than that on the non--

Concentration Field
When ξ u = 0.1, ξ l = 0.15, α = 3, β = 4, x = 3, λ = 1, κ = 2, ωt = 1, and Re ω = 0.1, Figure 9 shows the comparison of concentrations at the fixed point x = 3 of the microchannel with three different wall surfaces: non-roughness (C 0 ), in-phase (C + ), and out-of-phase (C − ).The concentration C + is higher than the concentration C 0 of the non-roughness, which indicates that the solute rate on the in-phase walls is higher than that on the nonroughness walls.On the contrary, the concentration C − is smaller than the concentration C 0 ; that is, the solute diffusivity on the out-of-phase walls is smaller than that of the nonroughness walls.It can be seen from Figure 8a,b that the velocity U + of the in-phase walls will increase with the increase in the roughness.However, the velocity U − of the out-of-phase walls will decrease accordingly.These lead to the difference in the solute diffusion of the three walls.Due to the shape of the out-of-phase walls, the fluid resistance is generated in the microchannel, which hinders the flow of the fluid [31].Consequently, solute diffusion in the microchannel is reduced.Figure 10 illustrates the changes in the concentration with the oscillation Reynolds number.It can be found that, with the increment of the oscillation Reynolds number, the local concentration C ± increases when the oscillation Reynolds number is below a certain critical value, and the concentration C ± will decrease when it exceeds this critical value.The reason is clear: although the large oscillation Reynolds number will improve the oscillation of the velocity, it will also lead to a decrease in the amplitude of the velocity.Figure 10 illustrates the changes in the concentration with the oscillation Reynolds number.It can be found that, with the increment of the oscillation Reynolds number, the local concentration C ± increases when the oscillation Reynolds number is below a certain critical value, and the concentration C ± will decrease when it exceeds this critical value.The reason is clear: although the large oscillation Reynolds number will improve the oscillation of the velocity, it will also lead to a decrease in the amplitude of the velocity.
Figure 11 shows the change in solute diffusions with the modulation charged parameters α and β.In Figure 11, α = β = 0 represents the microchannel with no modulated charged surface; i.e., only the effect of the symmetric walls potential is present.When α and β are not zero, they represent modulated charged surface microchannels.In these cases, the local concentration C + also increases due to the effect of the modulated charged surface.This is because modulating the charged surface increases the amplitude of the velocity, which in turn leads to faster diffusion of the solute.Due to the combined influence of roughness and modulated charged surface, the profile of concentration C − shows an oscillating nature.
Figure 12 illustrates the change in concentrations with the Péclet number Pe.In Figure 12, the diffusion of solute in the microchannel slows down as the Pe increases.This is because a low Péclet number Pe represents the diffusion term dominant.The flow with a small Reynolds number is considered in this study.Solutes in low Reynolds number flows are mainly affected by diffusion mechanisms.As the Péclet number Pe increases, the restricted diffusion mechanism causes a decrease in the local concentration C ± .
Figure 10 illustrates the changes in the concentration with the oscillation Reynolds number.It can be found that, with the increment of the oscillation Reynolds number, the local concentration C ± increases when the oscillation Reynolds number is below a certain value, and the concentration C ± will decrease when it exceeds this critical value.The reason is clear: although the large oscillation Reynolds number will improve the oscillation of the velocity, it will also lead to a decrease in the amplitude of the velocity.

Conclusions
In this paper, the AC EOF and mass transport in microchannels with modulated charged surfaces and roughness are investigated.The effects of modulating a charged

Conclusions
In this paper, the AC EOF and mass transport in microchannels with modulated charged surfaces and roughness are investigated.The effects of modulating a charged surface, sinusoidal roughness, and AC electric field on electrical potential, velocity, and mass transport are discussed in detail.The variations in potential, velocities, and concentration depend on some dimensionless parameters, such as oscillating Reynolds number (Re ω ), the roughness parameter (δ), the modulated charged surfaces parameter (α, β), and the Péclet number (Pe).By employing a theoretical analysis and a MATLAB R2021b calculation, the following conclusions are achieved: • When the modulated charge surface parameters (α, β) increases, the vertical velocity and the circulation flow are generated; the amplitudes of velocities also increase.In these cases, the local concentration of in-phase walls also increase.• As the roughness parameter δ amplifies, the oscillation of the vertical velocity is enhanced.The in-phase walls are conducive to the diffusion of solutes.• With the increment in the oscillation Reynolds number, the oscillations of the transverse velocity also increase.It can be found that the local concentrations will increase when the oscillation Reynolds number is below a certain critical value.

Appendix B. Numerical Algorithm
Utilizing the finite difference method solves the difference Equations ( 52)-(54).Since the origin of coordinates is taken at one end of the channel, the values of x and y are in the range of [0, 6] and [−1, 1].Consider x i (0 ≤ i ≤ n) and y j (0 ≤ j ≤ m) as discrete grid points.The steps τ= T k , h x = L N , and h y = 2H m divide the x and y axes into n and m components.The dimensionless concentration C of a point (x i , y j ) is denoted as C i,j .
The finite difference scheme is used in this paper as the second-order precision in space and the first-order precision in time; that is, O(τ + h2).The time step dt = 0.0005, the space step dx = dy = 0.1, and the Péclet number Pe varies from 0.3 to 1.The stability condition of the finite difference scheme at this time is as follows: The central difference scheme is as follows [40]:

Figure 1 .
Figure 1.Schematic diagram of mathematical model with roughness microchannel.

Figure 1 .
Figure 1.Schematic diagram of mathematical model with roughness microchannel.
Figure3shows the distribution of streamlines with the sinusoidal roughness amplitude δ, modulated parameters α, β, and Re ω under the conditions of λ = 1, κ = 2, ξ = 0.1, ξ = 0.15.Figure3adepicts streamlines in the non-rough microchannels with asymmetric wall potentials.The flow at this point is un-directional.Figure3b-ddescribe streamlines of the AC EOF for the in-phase walls, while Figure 3e-g depict streamlines of the AC EOF for the out-of-phase walls.Physically, the roughness and modulated charged surfaces change the ion distribution in the electrolyte solution in the microchannel, which affects the distribution of the potential in the EDL.Therefore, it will generate an influence on the flow of the

Figure 7
Figure 7 depicts the change in the transient dimensionless velocity  ± .In Figure7a,c, Re ω = 0.1, and in Figure7b,d, Re ω = 70.With the increase in the oscillatory Reynolds number Re ω , the velocity distribution gradually exhibits oscillatory characteristics.The velocity  ± oscillates in the region close to the wall.And it tends to zero in the region far from the wall at the large oscillatory Reynolds number.This is because, at high electric field frequencies, the diffusion time is greater than the oscillation period.

Figure 7
Figure 7 depicts the change in the transient dimensionless velocity  ± .In Figure7a,c, Re ω = 0.1, and in Figure7b,d, Re ω = 70.With the increase in the oscillatory Reynolds number Re ω , the velocity distribution gradually exhibits oscillatory characteristics.The velocity  ± oscillates in the region close to the wall.And it tends to zero in the region far from the wall at the large oscillatory Reynolds number.This is because, at high electric field frequencies, the diffusion time is greater than the oscillation period.

Figure 7
Figure 7 depicts the change in the transient dimensionless velocity U ± .In Figure7a,c, Re ω = 0.1, and in Figure7b,d, Re ω = 70.With the increase in the oscillatory Reynolds number Re ω , the velocity distribution gradually exhibits oscillatory characteristics.The velocity U ± oscillates in the region close to the wall.And it tends to zero in the region far from the wall at the large oscillatory Reynolds number.This is because, at high electric field frequencies, the diffusion time is greater than the oscillation period.Figure8illustrates the effect of the charge modulation surface and the sinusoidal roughness on the velocities U ± and V ± .In Figure8a,b, the dimensionless velocity U ± profile is presented, while Figure8c,d displays the dimensionless velocity V ± profile.Figure8a,bindicates that the amplitude of velocity U ± becomes enhanced with the application of modulation potential.This is because as the amplitude of the modulation potential heightens, the potential in the EDL also increases, leading to a rise in the amplitude of the velocity U ± .When the roughness increases, it can be observed from Figure8athat the velocity U + becomes incremental.However, Figure8bshows that the velocity U − slows down with the expansion of the roughness.In Figure8c,d, it is observed that as the amplitude of the modulation potential increases, the amplitude of the velocity V ± also increases.The rise of the roughness highlights the oscillation characteristics of the velocity V ± .

Figure 12
Figure 12 illustrates the change in concentrations with the Péclet number Pe.In Figure 12, the diffusion of solute in the microchannel slows down as the Pe increases.This is because a low Péclet number Pe represents the diffusion term dominant.The flow with a small Reynolds number is considered in this study.Solutes in low Reynolds number flows are mainly affected by diffusion mechanisms.As the Péclet number Pe increases, the restricted diffusion mechanism causes a decrease in the local concentration C ± .

Figure 12
Figure 12 illustrates the change in concentrations with the Péclet number Pe.In Figure 12, the diffusion of solute in the microchannel slows down as the Pe increases.This is because a low Péclet number Pe represents the diffusion term dominant.The flow with a small Reynolds number is considered in this study.Solutes in low Reynolds number flows are mainly affected by diffusion mechanisms.As the Péclet number Pe increases, the restricted diffusion mechanism causes a decrease in the local concentration C ± .